Crystalline character of high-magnetic-field cusp states in quantum dots 
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Conditional probability distributions from exact diagonalization are used to investigate the crys- 
talline or liquid character of the downward cusp states in parabolic quantum dots (QD's) at high 
magnetic fields. These states are crystalline in character for fractional fillings covering both low 
and high values, unlike the liquid Jastrow-Laughlin wave functions, but in remarkable agreement 
with the rotating- Wigner-molecule ones [Phys. Rev. B 66, 115315 (2002)]. The cusp states are 
precursors to the bulk fractional quantum Hall states (and not to the bulk Wigner crystal), since 
the collective rotation stabilizes the rotating Wigner molecule (formed in the QD) relative to the 
static one. 

PACS numbers: 73.21. La; 71.45.Gm; 71.45.Lr 
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The excitation energy spectra of two-dimensional N- 
electron semiconductor quantum dots (QD's) in high 
magnetic fields (B) exhibit downward cusps 0,[3>0>0 at 
certain magic angular momenta, corresponding to states 
with enhanced stability. These cusp states have been 
long recognized 0. OJ. Q as the finite-size precursors 
of fractional quantum Hall states in the bulk. In the 
literature of the fractional quantum Hall effect (FQHE) , 
ever since the celebrated paper |jj by Laughlin in 1983, 
the cusp states have been considered to be the antithesis 
of the Wigner crystal and to be described accurately by 
liquid-like wave functions, such as the Jastrow-Laughlin 
[111 (JL) and composite-fermion HG3 ( CF ) ones - This 
view, however, has been recently challenged by the 
explicit derivation of trial wave functions for the cusp 
states that are associated with a rotating Wigner (or elec- 
tron) molecule, RWM. As we discussed |4j earlier, the 
RWM wave functions, which are by construction crys- 
talline in character, promise to provide a simpler, but yet 
improved and more consistent description of the proper- 
ties of the cusp states, in particular for high angular mo- 
menta (L) [corresponding to low fractional fillings, since 
v = N(N — 1)/(2L)]. 

Issues pertaining to the liquid or crystalline character 
of the cusp states are significant in both the fields of QD's 
and the FQHE. Since the many-body wave functions in 
the lowest Landau level (high B) obtained from exact 
diagonalization (EXD), the RWM, and the CF/JL have 
good angular momenta L = Lq = N(N — l)/2 6], their 
electron densities are circularly symmetric. Therefore in- 
vestigation of the crystalline or liquid character of these 
states requires examination of the conditional probability 
distributions (CPD's, i.e., the fully anisotropic pair cor- 
relation functions). These calculations were performed 
here under high magnetic field conditions for QD's (in a 
disk geometry [llj) with N = 6 — 9 electrons, and for 
an extensive range of angular momenta. This allowed 
us to conclude that in all instances examined here (cor- 
rresponding to 0.467 > v > 0.111) the cusp states ex- 
hibit an unmistakably crystalline character, in contrast 



to the long held perception in the FQHE literature, with 
the RWM yielding superior agreement with the exact- 
diagonalization results. Furthermore, the RWM states 
are found to be energetically stabilized (i.e., exhibit gain 
in correlation energy) with respect to the correspond- 
ing static (symmetry-broken) Wigner molecules, from 
which the multideterminantal RWM wavefunctions are 
obtained through an angular-momentum projection |l2| . 
The CPD's are defined by 



P(r,ro) = ($|^£( ri 



rM ri -r )|$)/($|$), (1) 



where $(ri, ra, rjy) denotes the many-body wave func- 
tion under consideration. P(r,ro) is proportional to the 
conditional probability of finding an electron at r under 
the condition that a second electron is located at rg. This 
quantity positions the observer on the moving (intrinsic) 
frame of reference specified by the collective rotations 
that are associated with the good angular momenta of 
the cusp states. 

The CPD's for cusp states corresponding to a lower 
filling factor than v — 1/5, calculated for N = 6 electrons 
with L = 135 [y = 1/9, left column) and for N = 6 with 
L = 105 (u = 1/7, right column) are displayed in Fig. 1. 
Fig. 2 displays the CPD's for the cusp states with N = 6 
electrons and L = 75 (y — 1/5, left column) and N = 7 
and L = 105 [y — 1/5, right column). In both figures, the 
top row depicts the RWM case. The EXD case is given 
by the middle row, while the CF cases (which reduce to 
the JL wave functions for these fractions) are given by 
the bottom row. 

There are three principal conclusions that can be 
drawn from an inspection of Figs. 1 and 2: 

(I) The character of the EXD states is unmistakably 
crystalline with the EXD CPD's exhibiting a well devel- 
oped molecular polygonal configuration [(1,5) for N = 6 
and (1,6) for N = 7, with one electron at the center], in 
agreement with the explicitly crystalline RWM case. 

(II) For all the examined instances (covering the frac- 
tional fillings 1/9, 1/7, and 1/5), the JL wave func- 
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FIG. 1: Conditional probability distributions at high B for 
TV = 6 electrons with L — 135 (y = 1/9, left column) and 
L = 105 (y = 1/7, right column). Top row: RWM case. 
Middle row: The case of exact diagonalization. Bottom row: 
The Jastrow-Laughlin case. It is apparent that the exact di- 
agonalization and RWM wave functions have a pronouned 
crystalline character, corresponding to the (1,5) polygonal 
configuration of the rotating Wigner molecule. In contrast, 
the Jastrow-Laughlin wave functions fail to capture this crys- 
talline character, exhibiting a rather "liquid" character. The 
observation point (identified by the missing electron hump) 
was placed at the maximum of the outer ring of the radial 
electron density Q of the EXD wave function, namely at 
r = 7.318Z S for L = 135 and r = 6.442/ s for L = 105. 
The EXD Coulomb interaction energies (lowest Landau level) 
are 1.6305 and 1.8533 e 2 / kIb for L = 135 and L = 105, re- 
spectively. Here, Is = (hc/eB) 1 ^ 2 . The errors relative to 
the corresponding EXD energies and the overlaps of the trial 
functions with the EXD ones are: (I) For L = 135, RWM: 
0.34%, 0.860; JL: 0.50%, 0.665. (II) For L = 105, RWM: 
0.48%, 0.850; JL: 0.46%, 0.710. 



tions fail to capture the intrinsic crystallinity of the EXD 
states. In contrast, they represent "liquid" states in 
agreement with an analysis that goes back to the original 
papers @, by Laughlin. In particular, Ref. investi- 
gated the character of the JL states through the use of a 
pair correlation function [usually denoted by g(R)] that 
determines the probability of finding another electron at 
the absolute relative distance R = |r — ro| from the ob- 
servation point ro. Our anisotropic CPD of Eq. (JTJ is 
of course more general (and more difficult to calculate) 
than the g(R) function of Ref. 8]. However, both our 
P(r, r ) (for N — 6 and TV = 7 electrons) and the g(R) 
(for TV = 1000 electrons, and for v = 1/3 and v = 1/5) 
in Ref. reveal a similar characteristic liquid-like short- 
range-order behavior for the JL states, eloquently de- 
scribed in Ref. @ (see p. 249 and p. 251). Indeed, we 
remark that only the first-neighbor electrons on the outer 
rings can be distinguished as separate localized electrons 
in our CPD plots of the JL functions [see Fig. 1 and Fig. 



2]- 

(Ill) The pronounced crystallinity of the EXD states 
at the rather high v = 1/5 value (see Fig. 2) does not 
support the long held expectation 

'0S El (based on 
extrapolations of the JL wave functions to the bulk and 
comparisons with the static bulk Wigner crystal) that a 
liquid-to-crystal phase transition may take place only at 
the lower fillings with v < 1/7. 

Following the conclusion that the crystalline charac- 
ter of the cusp states in QD's is already well devel- 
oped for fractional fillings with the unexpected high value 
of v = 1/5, a natural question arises concerning the 
presence or absence of crystallinity in cusp states cor- 
responding to higher fractional fillings, i.e., states with 
1/5 < v < 1/3, and even with 1/3 < v < 1. To answer 
this question, we show in Fig. 3 the CPD's for the RWM 
(left column) and EXD (right column) wave functions for 
the case of TV = 7 and L = 45 (1/3 < v = 7/15 < 1), 
of TV = 8 and L = 91 (1/5 < v = 4/13 < 1/3), and of 
N = 9 and L = 101 ( 1/3 < v = 36/101 < 1). Unlike the 

N=6, L=75 (v= 1/5) N=7, L=105 (v= 1/5) 
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FIG. 2: Conditional probability distributions at high B for 
TV = 6 electrons and L = 75 [y = 1/5, left column) and 
for TV = 7 electrons and L = 105 (again v — 1/5, right 
column). Top row: RWM case. Middle row: The case of 
exact diagonalization. Bottom row: The Jastrow-Laughlin 
case. The exact diagonalization and RWM wave functions 
have a pronouned crystalline character, corresponding to the 
(1,5) polygonal configuration of the RWM for TV = 6, and 
to the (1,6) polygonal configuration for TV = 7. In contrast, 
the Jastrow-Laughlin wave functions exhibit a characteris- 
tic liquid profile that depends smoothly on the number TV of 
electrons. The observation point is located at ro = 5.4311b 
for TV = 6 and L = 75 and r = 5.883/ B for TV = 7 and 
L — 105. The EXD Coulomb interaction energies (lowest 
Landau level) are 2.2018 and 2.9144 e 2 / ' nl B for TV = 6, L = 75 
and TV = 7, L — 105, respectively. The errors relative to the 
corresponding EXD energies and the overlaps of the trial func- 
tions with the EXD ones are: (I) For TV = 6, L = 75, RWM: 
0.85%, 0.817; JL: 0.32%, 0.837. (II) For TV = 7, L = 105, 
RWM: 0.59%, 0.842; JL: 0.55%, 0.754. 
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FIG. 3: Additional CPD's at high B. RWM results: Left 
column. Results from exact diagonalization are depicted on 
the right column. Top row: N = 7 electrons and L = 45 
(1/3 < v = 7/15 = 0.467 < 1) Middle row: N = 8 electrons 
and L = 91 (1/5 < v = 4/13 = 0.308 < 1/3). Two bottom 
rows: N = 9 electrons and L = 101 (1/3 < v = 36/101 = 
0.356 < 1, see text for explanation). Even for these low magic 
angular momenta (high fractional fillings), both the exact- 
diagonalization and RWM wave functions have a pronouncd 
crystalline character [corresponding to the (1,6), (1,7), and 
(2,7) polygonal configuration of the RWM for N — 7, 8, and 
9 electrons]. The observation point is located at ro — 3.7761b 
for N = 7,L = 45, r = 5.W5l B for N = 8,L = 91, and r = 
5.2181b (outer) and r = 1.6621b (inner) for N = 9, L = 101. 



long held perceptions in the FQHE literature (which were 
reasserted in a recent manuscript the CPD's in Fig. 
3 demonstrate that the character of the cusp states with 
high fractional fillings is not necessarily "liquid-like" . In- 
stead, these high-i/ cusp states can exhibit a well devel- 
oped crystallinity associated with the (1,6), (1,7), and 
(2,7) polygonal configurations of the RWM, appropriate 
for N — 7, 8, and 9 electrons, respectively. Similar results 
were also found for the case of N = 6 electrons. Note 
that the case of N = 9 represents the smallest number of 
electrons with a non-trivial concentric-ring arrangement, 
i.e., the inner ring has more than one electrons. As the 
two CPD's [reflecting the choice of taking the observa- 
tion point (i"o in Eq. l[T]l) on the outer or the inner ring] 
for N = 9 reveal, the polygonal electron rings rotate in- 
dependently of each other. Thus, e.g., to an observer 
located on the inner ring, the outer ring will appear as 
uniform. 



Our two-step method for deriving the RWM wave func- 
tions is anchored in the distinction [15J between a static 
and a rotating Wigner molecule, with the collective rota- 
tion stabilizing the latter relative to the former. Further 
elaboration on this point requires generation of global 
ground states out of the cusp states, achieved through 
inclusion 0, Ha) of an external parabolic confinement 
(of frequency ojq). In the two-step method, the static 
WM is first described by an unrestricted Hartree-Fock 
(UHF) determinant that violates the circular symmetry 
[la ] . Subsequently, the collective rotation of the WM is 
described by a post-Hartree-Fock step of restoration of 
the broken circular symmetry via projection techniques 
[iH. We note that, in the limit A ^ oo, the static WM 
of the UHF develops to the bulk Wigner crystal 01 an d 
its more sophisticated variants fl3j ]. 

In general, the localized broken symmetry orbitals of 
the HF determinant are determined numerically via a 
selfconsistent solution of the UHF equations ^(|. Since 
we focus here on the case of high B, we can approximate 
the UHF orbitals (first step of our procedure) by (pa- 
rameter free) displaced Gaussian functions; namely, for 
an electron localized at Zj, we use the orbital 



with z = x + iy, Zj = Xj + iYj, and A = y/h/m*fl; 
SI = y/ojQ + where ui c — eB/{m*c) is the cyclotron 

frequency. The phase guarantees gauge invariance in the 
presence of a perpendicular magnetic field and is given in 
the symmetric gauge by ip(z, Zj; B) — (xYj — yXj)/2l 2 B , 
with Ib — y/hc/eB. We only consider the case of fully 
polarized electrons, which is appropriate at high B. 

We take the Zj's to coincide with the equilibrium po- 
sitions (forming nested regular polygons) of classical 
point charges inside an external parabolic confinement of 
frequency luq, and proceed to construct the UHF determi- 
nant ^ UHF [z] out of the orbitals u(^,Zj)'s, i = 1, ...,N. 
Correlated many-body states with good total angular 
momenta L can be extracted ^3] from the UHF determi- 
nant using projection operators. We stress that while the 
initial trial wave function of the UHF equations consists 
of a single determinant, the projected wave function is 
a linear superposition of many determinants. The pro- 
jected energy, corresponding to the symmetry-restored 
state with angular momentum L, is given by [l2j 

Epkj(L) = J 2 J M 7 )e^ L d 7 / fj «(7)e l7i d7, (3) 

with h(j) = (* UHF (0)|#|* UHF (7)) and rc( 7 ) = 
(* UHF (0)|* UHF ( 7 )), where f UHF ( 7 ) is the original UHF 
determinant rotated by an azimuthal angle 7 and H is 
the many body Hamiltonian (comprised |l2| of the ki- 
netic energy with a vector potential for the magnetic 
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FIG. 4: Stabilization energies A£| s ain = £§ S RJ - £?uhf for 
TV = 6 (dashed curve) and N = 7 (solid curve) fully polarized 
electrons in a parabolic QD as a function of B. The troughs 
associated with the major fractional fillings (1/3, 1/5, and 
1/7) and the corresponding ground-state angular momenta 
are indicated with arrows. We have extended the calculations 
up to B = 120 T (not shown), and verified that A£f s ain re- 
mains negative while its absolute value vanishes as B — » oo. 
The choice of parameters is: hu>o = 3 meV (parabolic confin- 
ment), m* — 0.067m e (electron effective mass), and k = 12.9 
(dielectric constant). 



field, an external confinement, the Coulomb interelectron 
repulsion, and the Zeeman term). The UHF energies are 
simply given by -Etjhf = h(0)/n(0). 

We note that, unlike the UHF ground state (describing 
a static Wigner molecule) which has no good angular mo- 
mentum, the ground states of the RWM exhibit good an- 
gular momenta (labeled as L gs ) that coincide with magic 
ones [we denote the ground-state energy of the RWM as 

The stabilization of the rotating WM relative to the 
static one [namely the fact that -Eppj < -E-uhf, see Fig. 
4] is a purely quantum effect. This energy gain demon- 
strated here for N — 6 and 7 electrons, is in fact a gen- 
eral property of states projected out of trial functions 
with broken symmetry. This is due to a theorem [lsf 
stating that at least one of the projected states (i.e., the 
ground state) has an energy lower than that of the orig- 
inal broken-symmetry trial function (e.g., the UHF de- 
terminant considered above). 

While our focus here is on the behavior of trial and 
exact wave functions in (finite) QD's in high magnetic 
fields, it is natural to inquire about the implications of 
our findings to FQHE systems in the thermodynamic 
limit. We recall that appropriate trial wave functions for 
clean FQHE systems must possess a good angular mo- 
mentum L > Lq, a property shared by both the CF/JL 
and RWM functions. We also recall the previous find- 
ing that for small fractional fillings v the static Wigner 
crystal (which is a broken-symmetry state with no g ood 
angular momentum) is energetically favored pi y| [lfj, Il3| 
compared to the (liquid-like) CF/JL wave function. In 



contrast, the RWM wave functions remain lower in en- 
ergy than the corresponding static crystalline state for 
all values of N and v, even in the thermodynamic limit. 
This is due to the fact that the aforementioned "energy- 
gain" theorem 18] applies for any number of electrons N 
and for all values of the magnetic field B [l9( . 

Since the rotating Wigner crystal carries a current, 
while the static crystal is insulating, we may conjecture 
that a transition at lower fractional fillings from a FQHE 
conducting state to an insulating Wigner crystal can- 
not occur spontaneously for clean systems. Therefore, it 
should be possible to observe (see, e.g., Ref. |2(J) FQHE 
behavior at all fractional fillings in a clean system. In 
practice, however, impurities and defects, may pin the 
RWM, so that one of the main challenges for FQHE ob- 
servation at such low fillings relates to fabrication of high 
mobility samples 0, • 

In summary, we have carried out the first systematic 
investigations (for 6 < N < 9) of cusp states in parabolic 
quantum dots at high magnetic fields. Our anisotropic 
conditional probability distributions from exact diago- 
nalization show that these states are crystalline in char- 
acter for both low and high fractional fillings, unlike 
the liquid-like Jastrow-Laughlin 0, @ wave functions, 
but in remarkable agreement with the recently proposed 
rotating- Wigner-molecule Q ones. The cusp states of 
finite- N parabolic QD's are precursors to the bulk frac- 
tional quantum Hall states (and not to the bulk Wigner 
crystal) due to the collective rotation, which stabilizes 
the rotating Wigner molecule (having a good angular mo- 
mentum) relative to the static one (that exhibits broken 
symmetry) . 

We thank M. Pustilnik for comments on the 
manuscript. This research is supported by the U.S. 
D.O.E. (Grant No. FG05-86ER45234). 



[1] R.B. Laughlin, Phys. Rev. B 27, 3383 (1983). 
[2] P.A. Maksym and T. Chakraborty, Phys. Rev. Lett. 65, 
108 (1990). 

[3] J.K. Jain and T. Kawamura, Europhys. Lett. 29, 321 
(1995). 

[4] C. Yannouleas and U. Landman, Phys. Rev. B 68, 035326 

(2003); 66, 115315 (2002). 
[5] S.R.E. Yang, A.H. MacDonald, and M.D. Johnson, Phys. 

Rev. Lett. 71, 3194 (1993). 
[6] A.L. Jacak, P. Hawrylak, and A. Wojs, Quantum Dots 

(Springer, Berlin, 1998), in particular Ch. 4.5. 
[7] R.B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 
[8] R.B. Laughlin, in The Quantum Hall Effect, Edited by 

R.E. Prange and S.M. Girvin (Springer, New York, 1987) 

p. 233 

[9] J.K. Jain, Phys. Rev. B 41, 7653 (1990). 
[10] S.S. Mandal, M.R. Peterson, and J.K. Jain, Phys. Rev. 

Lett. 90, 106403 (2003). 
[11] The results in Figs. 1-3 were obtained without consider- 



5 



ation of a confinement, and they are not influenced (high 
B) by inclusion of a parabolic confinement plPlll For 
the influence of the confinement on the energies, see Refs. 

Sana. 

[12] C. Yannouleas and U. Landman, J. Phys.: Condens. Mat- 
ter 14, L591 (2002); Phys. Rev. B 68, 035325 (2003), and 
references therein. 

[13] P.K. Lam and S.M. Girvin, Phys. Rev. B 30, 473 (1984). 

[141 G.S. Jeon, C. -C. Chang, and J.K. Jain, 
cond-mat/0310287 

[15] C. Yannouleas and U. Landman, cond-mat/0311480 

[16] C. Yannouleas and U. Landman, Phys. Rev. Lett. 82, 
5325 (1999); Phys. Rev. B 61, 15 895 (2000). 

[17] D.J. Yoshioka and P.A. Lee, Phys. Rev. B 27, 4986 
(1983). 



[18] See sect. 3 in P.-O. Lowdin, Rev. Mod. Phys. 34, 520 
(1962). 

[19] We remark that, at present, computational difficulties 
prohibit a straightforward extrapolation of energy to the 
thermodynamic limit as a function of N and v. This is 
the case especially for the disk geometry used here, which 
is advantageous for studies of crystalline structures, com- 
pared to the computationally easier spherical geometry. 

[20] W. Pan et al, Phys. Rev. Lett. 88, 176802 (2002). 

[21] L. Pfeiffer and K.W. West, Physica E 20, 57 (2003). 

[22] We remark, however, that the stabilization energy dimin- 
ishes as -B — > oo (see Fig. 4), and as a result the impurities 
become more efficient in pinning the WM for the lower 
fractional fillings (i.e., higher angular momenta). 



